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Abstract 

Description of a nuclear system in its ground state and at low exci¬ 
tations based on the equation of state (EoS) around normal density is 
presented. In the expansion of the EoS around the saturation point ad¬ 
ditional spin polarization terms are taken into account. These terms, 
together with the standard symmetry term, are responsible for appear¬ 
ance of the a-like clusters in the ground state conhgurations of the N=Z 
even-even nuclei. At the nuclear surface these clusters can be identified 
as alpha particles. A correction for the surface effects is introduced for 
atomic nuclei. Taking into account an additional interaction between 
clusters the binding energies and sizes of the considered nuclei are very 
accurately described. The limits of the EoS parameters are established 
from the properties of the a, ^He and t particles. 


1 Introduction 

In 1936 Bethe and Bacher [Ij and in 1938 Hafstad and Teller predicted 
that alpha particle structures 1^0 could be present in atomic nuclei. Convinc¬ 
ing arguments for existence of such structures were provided by systematics of 
the binding energies of the even-even nuclei with equal number of protons and 
neutrons (see also [^), as well as by systematics of the binding energy of the 
additional neutron in nuclei like ®B, and The latter systematics could 
be explained by assuming that the valence neutron moves in a multi-center 
potential with centers identified with alpha particles. Today, in a commonly 
accepted point of view one considers that alpha particle clusters can appear in 
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the nuclear matter at subsaturation densities. Such conditions are met e.g. at 
the nuclear surface where, in not too heavy nuclei, almost the whole nuclear 
mass is concentrated. 

Subtle changes of the ground state energy of atomic nuclei as a function of 
nucleon number can be described by the shell model, which neglects correlations 
between nucleons. The main idea of the present approach is to verify whether 
these subtle changes, in case of n-alpha nuclei, could be attributed to the ap¬ 
pearance of alpha-like clusters which modify the structure of the wave function 
and introduce strong correlations between nucleons. In such an approach the 
system in its ground state behaves like a crystal, with stationary configuration 
and shape and with definite distances between the wave function centers for a 
given nucleus. The subtle changes of the binding energy due to appearance of 
alpha-like clusters in the wave functions of n-alpha nuclei will be determined 
within an extended version of the model proposed in j^. 

It will be shown in the following that the hypothesis of alpha structures 
in the n-alpha nuclei can indeed describe the binding energy systematics of 
Bethe mi It will be also shown that the hypothesis of the alpha structures 
helps to understand the evolution of the root mean square (rms) radii and of 
the density profiles as a function of the number of alpha particles in the n-alpha 
nuclei. 

The paper is organized as follows: section 2 describes the adopted interac¬ 
tions in the nuclear matter. Section 3 presents constraints on the introduced 
parametrization resulting from the properties of light charged particles. Section 
4 is devoted to the ground state properties of n-alpha nuclei. First, it inves¬ 
tigates the effects of clusterization on the ground state configurations. Then, 
some corrections to the Hamiltonian, related to the many body character of the 
nuclear interaction and to the bosonic nature of the wave function are intro¬ 
duced. It will be shown that the introduced corrections make the structural 
corrections stronger. Finally, the model predictions for binding energies, rms 
radii and density profiles are presented and compared with the experimental 
values. Section 5 contains conclusions. 


2 Interactions in nuclear matter around satura¬ 
tion 

In order to predict correctly a state and evolution of a nuclear system a reliable 
description of the interaction between nucleons is needed. In case of a strong 
interaction it is not an easy task. As we know the fundamental theory governing 
the nucleon-nucleon interactions is the quantum chromodynamics (QCD), where 
interactions are calculated from the physics of quarks and gluons. Unfortunately, 
very limited progress has been achieved in application of QCD for the descrip¬ 
tion of interaction between quark-gluon clusters. Therefore a phenomenological 
approach, associated to some specific conditions has been applied. Since we are 
interested in the low-energy region in which nucleons are not excited internally. 
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one can treat nucleons as quasi elementary particles. Additionally, if we apply a 
non-relativistic approximation, the mutual nucleon interaction can be described 
by a respective potential. As can be found in the text books, details of the 
form of such a potential are deduced from nuclear interaction symmetries and 
mechanism of exchange of bosons. 

For a final test of this type of a potential and the determination of its pa¬ 
rameters one uses a comparison with the experimental data. Here the nucleon- 
nucleon scattering data or light nuclei properties are usually used. A general 
conclusion from this kind of research is that a simple two body potential is not 
sufficient to describe the behavior of the system of more than two nucleons. One 
should add at least the three-body force. In particular, in the nuclear system, 
the three-body force is acceptable because the strong interaction is generated 
by the exchange process of gluons. The n-body interaction can be interpreted 
as coming from the multi-pion exchange processes. 

These considerations show the difficulties in the description of even the sim¬ 
plest nuclear system. The extension of the description of interaction by the 
three-body forces shows that the exchanged particles (bosons) play a significant 
role in the interaction process. 

In order to describe dynamics of the system one has to express the change 
of energy associated with its evolution. In the description of energy through 
the potential one has to define its form. Hence the wave function of a system 
of identical particles must be either symmetric or antisymmetric and it must be 
taken into account when defining the energy of the ground state. In case of the 
fermionic systems one should use antisymmetric wave function which is a serious 
difficulty in the description of the collision dynamics. Some proposals to solve 
this problem are shown for instance in jTjj^. Other approaches, introducing 
some kinds of the momentum dependent Pauli potentials and Pauli-blocking 
schemes, which are meant to mimic the effects of antisymmetrization and are 
much less time consuming, are also proposed, see e.g. 
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Therefore, in the present approach, we propose a description in which parti¬ 
cles are distinguishable and energy associated with the fermionic motion affects 
the resulting energy density which is parametrized and used instead of the po¬ 
tential. Such a concept is present in the description of the liquid drop model 
(LDM) and in some density functional models, e.g. in [^. It is obvious that 
such an approach requires the knowledge of the local density of matter and 
the knowledge of the wave function describing the system. In the present de¬ 
scription we use a standard one, commonly applied in such cases, i.e. the wave 
function is defined as a simple product of minimal wave packets (eq. (1)), as 
e.g. in |16| . The main argument for disregarding the antisymmetrization in the 
present approach is the form of the proposed EoS which includes the spin and 
isospin terms. These terms force the unlike nucleons to form “bosonic” alpha¬ 
like clusters for which, presumably, the antisymmetrization may no longer play 
a decisive role. 

Each minimal wave packet in the adopted wave function has a defined spin 
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and a specific charge for a selected quantization axis: 




( 1 ) 



( 2 ) 


where cr'j.{r), (r^), (pfc), are the width (the position variance of the fc-th nucleon) 
of a Gaussian wave packet and the mean position and mean momentum of 
each of the A nucleons, respectively. Every partial wave function is labeled by 
Ik = n or Ik = p and Sk =t or Sk =)., denoting the isospin and spin (or 
precisely, their projections on a specified quantization axis) associated with a 
given nucleon. Variables (r^), (p^) and cr^(r) are, in general, time-dependent 
parameters describing the wave functions. 

In the present approach the nuclear matter is treated as a four component 
fluid characterized by the respective densities: 

Pp-f - for protons with spin up, 

Ppl - for protons with spin down, 

Pn^ - for neutrons with spin up, 

Pni - for neutrons with spin down. 

Density distributions in the system are uniquely determined by parameters of 
the wave function (mean position, mean momentum and variance of the position 
of each nucleon). The description of the system evolution due to the time 
evolution of these parameters is, in general, defined by the Dirac-Frenkel time- 
dependent variational principle |17| . This approach describes the interaction in 
a self-consistent way. 

Assuming that nucleons are moving inside the nuclear system, the average 
energy of the system is a sum of the average potential energy and of kinetic en¬ 
ergy associated with the fermionic motion. This is the case when one disregards 
the energy of the ordered motion, which for the low excitations of the system is 
considered later on as a correction (see formula @). 

According to the above idea we now present a method of determining the 
average energy associated with each nucleon. For the description of this average 
energy we assume the existence of a scalar field e (pp^, Ppi, Pni), determined 
by the local densities pp^(r), pp^(r), pn^r), p„^(r). 

If, for every k — th nucleon, the corresponding wave packet ^(t^ikSk determines 
the probability Pfc(r) = of finding a nucleon at a given point r then 

it is assumed that the average energy e/j associated with this nucleon is defined 
by the mean value and variance of the field e (pp^, Pp^, Pnt i can be 

expressed as: 


e* = (£)fc + (e) 


(3) 


where: 
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(5) 


o-fe (e) = J Pk{r) {s{pp^,ppi,pn^,pni) - {e)kf d^Y 

and A is a parameter related to the surface energy, similarly as in |18| . 

The EoS of the nuclear matter is defined for an infinite system, assuming its 
homogeneity and isotropy. For such a matter, variance of the associated field 
£ (PpT> Ppi^ Pni) vanishes and the average value of the energy per nucleon is 
determined by the mean value of the field alone. The field e {pp^, Pp^, pn^, pn\) 
is in this case equivalent to the EoS of the nuclear matter. Thus, the functional 
form of the field e {pp^, Pn-\, Pn^) can be obtained from the EoS by applying 

the local density approximation. 

In the present work we use a 12-parameter cubic form of the EoS proposed 
in . In the vicinity of the saturation density it can be approximated using a 
second order Taylor expansion: 


e — Coo + —C + 


where: 



(6) 

O 

I £ 

II 

(7) 

g _ Pn Pp 

P 

(8) 

Pnt “ 

Vn = 

P 

(9) 

Pp\ - Ppl 

'Ip — 

P 

(10) 


in which p and pg are the total nuclear matter density and the density of isospin 
and spin balanced matter at saturation, respectively. 

The first two terms in ^ are just a standard, 6-parameter form of the EoS 
approximation commonly used in the present day experimental and theoretical 
studies, see e.g. [^. The first term describes the symmetric matter in a bal¬ 
anced system (zero isospin and spin), with Kq being the compressibility of the 
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symmetric nuclear matter. The second one is a standard approximation of the 
symmetry energy: 


e/ = e/o + (11) 

with e/o being the Wigner constant and the coefficients L/ and Kj being the 
slope and curvature, respectively. 

The last two terms of ^ are the main novelty of the present approach (see 
also [^). In analogy to the isospin symmetry energy, they describe the spin 
symmetry energies for neutrons and protons separately: 

/ pr 

CiiO “I “I C (12) 

and the energy of the mutual, spin interaction of protons and neutrons: 

= (13) 

The terms with indices ii, i.e. the enQ, La and Ka, (constant, slope and 
curvature) describe the energy of protons or neutrons. Indices ij indicate that 
the symmetry energy refers to the mutual interaction of protons and neutrons. 
Determination of the expansion parameters in ^ will be discussed in the next 
chapters. 

Using the field e (pp^, Ppi, Pn^, Pni) one can express the average value of the 
Hamiltonian for a given system of A nucleons as: 




k—A / \2 

(P/e) 


k=A 

+ A <Jk{£) 


k=l 


($|Uc|^) (14) 


k=A 

k=l 

where tojv is the nucleon mass and the last term describes the Coulomb energy. 
As later on we will be searching for the ground state configurations of the 
even-even nuclei with N = Z, it can be assumed that the average momenta 
of the wave packets (p/,) are equal to zero and, thus, the first term in the 
expression (14) can be neglected. Following a standard approach (see e.g. [T^) 


it is assumed that the parametrization of the EoS includes implicitly the kinetic 
contribution coming from the Fermi motion. As has been demonstrated in |18| , 
the adopted cubic approximation of the EoS is flexible enough to reproduce 
typical predictions of sophisticated models of the nuclear matter which explicitly 
take into account the kinetic contribution resulting from the Fermi motion. 

Since the local density can be written as: 


k=A 


, (r) = ^ Pk (r), 


(15) 


the average value of the Hamiltonian can be expressed as: 

k=A 


($|iJ|$)= f e{pp^,ppi,pnt,Pni) Pir)cfr +X^akie) + {^\Vc\<^) (16) 

k=l 
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As one can see the Hamiltonian is a sum of three components, which can be 
interpreted as a volume, surface and Coulomb energies. Therefore the present 
description can be interpreted as a microscopic realization of the LDM. 

Finding the ground state configuration of a group of A nucleons within the 
adopted parametrization of the wave function and for given parameters of the 
EoS, consists in finding the values of parameters cr^(r), (r^), and (pfc), which 
minimize the value of Hamiltonian. The procedure of finding the ground state 
configurations has been described in Sect. 4 of [^. We would like to emphasize 
that the proposed parametrization of the EoS, with two more terms in (6), 
is responsible for driving the nucleons of different types, pf, pi, nf, ni, to 
form quadruplets, i.e. to group into a-like clusters during the minimization 
procedure. These two terms play a similar role as the standard symmetry term 
which attains a minimum when the neutron-proton asymmetry S of (7) tends to 
zero, i.e. when the neutron-proton couples group together. This a-like clustering 
property of (6), is demonstrated in Figure]^ 

The histogram in Fig. presents the evolution of the binding energy of 
12 nucleons during the process of searching for the ground state configuration 
which minimizes the model Hamiltonian (16). The 6 panels of the inset show 
“snapshots” of the projections of the mean positions of the nucleon wave packets 
onto the x-y plane, at selected stages of the minimization procedure. Panel 1 
presents the initial random configuration, panels 2-5 some intermediate ones 
and panel 6 shows the final ground state configuration of As can be seen, 
the minimum of the Hamiltonian is achieved for a configuration with three a- 
like clusters. These clusters appear in the corners of an equilateral triangle. 
Note, that in the figure the triangle is not perfectly equilateral due to the 2- 
dimensional projection. 

It has to be emphasized, that if clusters are allowed to appear in the con¬ 
sidered system the description has to be modified. Modifications of the Hamil¬ 
tonian due to clustering will be discussed in section 4.2. 


3 Constraints on the EoS from the properties of 
light charged particles 

Ground state properties of light charged particles, LCP, such as d, t, ^He or 
a can provide important constraints on parameters which describe interactions 
in the nuclear matter for at least four reasons. First, one can assume that in 
their structure clusters of matter (obviously beyond nucleons) can be formed. 
Therefore one may expect that clustering of matter has no effect on their binding 
energies and sizes. Second, these LCPs do not contain more than one a particle 
and thus one can neglect the a — a interactions (see Sect. 4.2) and use the 
simpler Liquid-Drop-like Hamiltonian (16). Third, the measured density profiles 
of these light particles are relatively well described by Gaussian distributions 
which is in line with the adopted parametrization of the wave function. Fourth, 
simplicity of these particles assures vanishing of some types of interactions in (6) 
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Figure 1: Evolution of the binding energy of 12 nucleons (Sxpf, SxpJ,, Sxnf and 
SxnJ,) during the minimization procedure. The insets 1-6 show the projections 
of the configurations of nucleons on the x-y plane at selected stages of the 
minimization. Panel 1 represents the random initial configuration. Panel 6 
shows the projection of the final ground state configuration of 


and selective increase of sensitivity to some other ones. For instance, parameters 
related to the isospin symmetry energy (second term in (|^) in principle do not 
affect the binding energy of a deuteron and alpha particle. Similarly, due to the 
saturation of the proton (neutron) spin interactions in ^He (t), the mutual spin 
interaction (fourth term in (|^) vanishes in their case. 

In case of LCPs we can assume that the average value of the Hamiltonian is 
completely determined by the field £ (pp^, Pnti Pnt)’ Coulomb inter¬ 

action and by the surface energy. This assumption may not be valid for more 
complex systems in which mutually interacting clusters can be formed, eg. for 
®Be—>■ a — a or ®Li —a — d. For such systems one can expect additional inter¬ 
actions between clusters due to the specific many-body forces. The mean value 
of the Hamiltonian (16) depends on 13 parameters: 12 parameters of the field 
e(r) plus the surface tension parameter A. Binding energies and rms radii of the 
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ground states for d, t, ^He and a provide 8 equations, thus insufficient number 
of constraints for a unique choice of parameters describing the interaction. In 
the present section we will try to deduce a subspace of this 13-parameter-space 
for which the above mentioned properties of LCPs are reproduced. 

3.1 Constraints on po, ego, Kg and A from the properties of 
an alpha particle 

It is clear that the greatest reduction of the number of parameters can be ex¬ 
pected in case of an alpha particle. Here, the spin and isospin dependent in¬ 
teractions are practically absent. Thus, properties of alpha particles should be 
described by four parameters only: pg, egg, Kg and A. Constraints on these 
parameters are imposed by the rms radius and the binding energy of an alpha. 
These parameters cannot be determined unambiguously. However, there is a 
flexibility in determining a two-dimensional hypersurface, in the space pg, egg, 
Kg and A, in which the searched values of these parameters exist. 

In order to determine this hypersurface the following variable is defined: 



where. Bag, Bam, fae, Tam the experimental and model values of the binding 
energy and rms radius, respectively. Note that a variable defined this way should 
be zero in the searched subspace. 

In order to determine the correlations among the parameters in a searched 
subspace, the variables pg and Kg are predefined on a grid, and the remaining 
two parameters are found by requiring that < 10~^ for the model ground 
state configurations of an alpha particle. The results are shown in Fig. 

The figure presents the constant saturation density curves as a function 
of the other three parameters specified on the axis, which ensure a correct 
binding energy and rms radius of an alpha particle for the model ground state 
configurations. The dotted and dashed lines show how one can determine the 
parameters A and egg for a given compressibility Kg and saturation density 
Pg. In the illustrated example the parameters pg — 0.185 fm“^ and Kg = 270 
MeV determine the values of A = 0.4 and eoo = —16.79 MeV. These values of 
the parameters have been suggested by additional studies of the competition 
between Coulomb, surface and volume energies in the binding energy calculated 
for the ground state configurations of the heavier nuclei (see Table 1). This 
problem will be considered later on. 

3.2 Constraints on eig, L/, Ki, eng, La and Kn from the 
properties of t and ^He 

After fixing the parameters pg, ego, Kg and A let us try to examine the other 
ones. For t and ^He the ground state energy does not depend on the choice of 
parameters eijg, Lij and Kij. This results from vanishing of the rjnVp product in 
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Figure 2: Constant saturation density cnrves, pQ, as a function of eoo, Kq and 
A, constrained by the ground state properties of an a particle. 


the last term of Thus, the binding energies and rms radii for t and ^He can 
provide four equations for six unknown EoS parameters: e/o, T/, Kj, eao, La 
and Ka, and so, two eqnations are still missing. The two free parameters can be 
constrained from the properties of heavier nuclei as it was done in the case of the 
fully symmetric nuclear matter. At this stage it is only possible to examine the 
snbspace defining the limits for the considered parameters from the properties 
of the t and ^He nnclei. We notice the fact that the parameters associated 
with isospin: e/o, Lj and Kj differentiate between the binding energy and the 
rms radii of t and ^He while the parameters eao-, La and Ka affect equally the 
energy and the size of these particles. 

Using a similar method as in the case of an alpha particle one can find 
relationships between the parameters e/o, Kj and L/, by requiring a precise 
description of the binding energies and radii of the t and ^He particles. Here, 
the values of e/o and Kj were varied independently on a grid and the value of 
L[ was found from the minimization of The resnlts of such calculations 

are shown in Fig. 

The figure presents a relationship between the parameters e/o, Lj and Kj. 
Selection of two of them allows to determine the third one. It shonld be stressed 
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Figure 3: Correlation between the slope, L/, and curvature, Ki, of the isospin 
symmetry energy for different Wigner energies, e/oi obtained from fitting the 
binding energies and rms radii of the t and ^He. 


that the above correlations have been obtained as a result of the procedure 
reproducing the t and ^He ground state properties with the parameters po = 
0.185 fm“^, Kq = 270 MeV, A = 0.4 and eoo = —16.79 MeV fixed from the 
ground state properties of an alpha particle (Fig. [^. Three other parameters: 
CizOj La and were found to be almost independent of e/o and Kj. 

3.3 Discussion of the parameter values 

The final values of the EoS parameters adopted in this work have been obtained 
by taking into account also the binding energy and size of a deuteron. The 
values have been summarized in Tables and For this example set of values 
the model reproduces exactly the experimental values of the rms radii and of 
the binding energies of the d, t, ^He and a particles (see also Figs. |^and[^. 

From the point of view of the present work, which focuses on the nuclear 
systems with even and equal numbers of protons and neutrons (i.e. the spin- 
balanced matter), the parameters: euo, La, Ku and e^Oj Lij^ Kij from the 
expansion ([^ are irrelevant. However, they do play an important role when 
constructing the ground state configurations and determining the properties of 
nuclei composed of asymmetric matter. 

As far as the interactions between neutrons and protons are concerned, the 
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Coulomb repulsion between protons causes that the proton matter tends to 
locate itself in the outer regions of the nucleus. This leads to differences between 
the neutron and proton density distributions, even for symmetric nuclei, and 
triggers additional interactions associated with the isospin symmetry energy 
(proton or neutron skins can be created). This effect can be negligible in case 
of light nuclei, while for heavier ones it may be more significant. Therefore, 
for heavier nuclei, one has to use a reliable set of isospin symmetry energy 
parameters e/o, T/ and Kj, even for symmetric nuclear matter. 

As has been shown above, the binding energies and the rms radii of t, ^He 
and a impose certain constraints on the parameters of the EoS. In case of 
the eoo, Po, Kq and A parameters, once two of them are fixed by some other 
constraints or means, or from some other evidences, the remaining ones can be 
determined using the obtained correlations (Fig. |^. Here, the compressibility 
parameter of symmetric matter Kq should fall within the range of about 250±50 


MeV according to the recent experimental constraints (flow interpretation 20 
monopole vibrations , subthreshold kaon production [Mpi] ). It has been 

fixed to 270 MeV in the present approximation. The volume energy parameter 
eoo should take the value of about -16 MeV according to the LDM. In the present 
approximation it has been assumed to have the value of -16.79 MeV. The values 
of the remaining two parameters, po A have been then uniquely determined 
from the correlation plot of Fig. 

In case of the symmetry energy, the choice of any two parameters out of e/o, 
L/, AT/, allows to determine the third one (Fig. [^. The compilation of the 
latest results of the nuclear physics experiments and astrophysical observations 
suggests that the value of the symmetry energy constant, e/o, and of the slope 
parameter, Lj, should fall within a range of about 32 ± 1 MeV and of about 
59 ± 17 MeV, respectively. In the present approximation the values of these 
parameters have been assumed to be 32 and 59.405 MeV, respectively. These 
values yield the value of Kj to be about -1250 MeV, according to the correlation 
plot of Fig. The experimental constraint of |26| on the symmetry energy 
compressibility parameter, Kj, yields its value to be about —50± 200 MeV. 
Another compilation locates the theoretical values of AT/ between —400 
and +466 MeV and its experimental values in a range from — 566± 1350 MeV 
to 34 ± 159 MeV. Thus, this parameter is least constrained experimentally 
and theoretically so far. Its value obtained within the model, preferring very 
soft symmetry energy, is consistent with the negative experimental constraint, 
thanks to its large error bar. However, we should stress that the adopted set 
of values is not at all unique. As seen from the correlation plots, some other 
choice of parameters is equally well possible. It should be also emphasized that 
the assumed value of po equal to 0.185 fm“^, which fits well properties of the 
LCPs, and is higher than the commonly used values of 0.16 — 0.17 fm“^, may 
influence values of the remaining parameters. 

In order to describe the exotic nuclei and provide more precise parameters 
of the symmetry energy, definitely some more advanced parameter searches will 
be needed, involving e.g. fitting the properties of the neutron skin nuclei and 
also of those approaching the proton drip line across the chart of the nuclides. 
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but this is beyond the scope of the present work. 


Po 

eoo 

Ko 

e/o 

Li 

Ki 

A 


1 

fm'^ 


[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

1 


0.185 

-16.79 

270 

32 

59.405 

-1250 

0.4 


Table 1: Parameters of the standard EoS which together with a parameter of 
the surface energy A are reproducing the ground state energy and rms radius of 
the alpha particle, t and ^He. 




Ku 

S-ijO 

Lij 


[MeV] 


1 

fm'-> 


[MeV] 

[MeV] 

[MeV] 

[MeV] 

82.37 

242.48 

19.909 

-1.5 

182. 

1210. 


Table 2: Parameters of the spin interaction chosen in such a way that together 
with parameters from Table 1 they reproduce binding energies and sizes of d, t 
and ^He. 

A concise summary of the procedure of fixing the parameter values is pre¬ 
sented in the Appendix. 


4 Energy and size of nuclei with equal and even 
number of neutrons and protons 

A method of finding the ground state configuration based on minimization of the 


average value of the liquid drop-like Hamiltonian (16) with respect to the wave 


function parameters (r^) and ak{r) with the (pfc) = 0 constraint, is described in 
H- The ground state configuration parameters determined this way indicate 


that the corresponding wave function has some numbers of centers, which can 
be interpreted as the centers of alpha-like clusters in the nuclear matter. Some 
features of these ground state configurations have been discussed in without 
taking into account interactions and correlations among clusters. In the present 
approach modification of the ground state configurations due to the presence of 
the alpha-alpha interactions will be investigated. 

The main aim is to identify such symmetries of the ground state configura¬ 
tion, and thus of the corresponding wave function, which could be responsible 
for the experimentally observed changes of the binding energies and rms radii of 
the n-alpha nuclei as a function of n. We would like to check whether it is pos¬ 
sible to identify specific configurations related to the magic numbers. In other 
words, we would like to verify whether the shell corrections necessary to describe 
deviations from the smooth predictions of the LDM around the magic numbers, 
could be interpreted in terms of a change of the ground state configuration re¬ 
sulting from the corrections to the Hamiltonian related to the cluster-cluster 
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interactions. As will be justified in the following, the ground state configura¬ 
tions corresponding to the local minima of the Hamiltonian can be classified 
according to the nuclear charge, Z, in the following way: 

i. Z< 8 The structures are identical with the structures proposed 

in and also in and are characterized by no alpha 
cluster in the center of the system. 

ii. 8 <Z< 20 A core represented by one alpha-like cluster appears in the 

center of the system. 

iii. 20 <Z< 28 A core with a mass corresponding to two alpha particles 

appears in the center. 

The structures for cases i and ii are visualized in Fig. |^(see also [^). 

1 < Z < 8 8 < Z < 20 


perspective view projections relative 

rotation 



Figure 4: Proposed structures of n-alpha nuclei. These structures ensure high 
level of symmetry for a given number of alpha particles, (see 0 )- The density 
distribution of the central cluster (the red one) has quite a big spread and 
therefore it cannot be identified as an alpha particle. 


The structures presented in Fig. [^result from the analysis of a large number 
of configurations obtained within a procedure of minimizing the Hamiltonian 
(16) when starting from a random configuration (see example in Fig. [^. 
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Majority of thus obtained configurations could be classified as those presented 
in the previous paper j^, with alpha-like clusters located mostly at the surface 
and with a moderate level of symmetry. Binding energies for these structures 
followed the smooth LDM trend and could not reproduce the “fine” structure 
of binding energies around the magic numbers. In some cases however, the 
minimization procedure did not yield the global minimum, but ended up in 
some local one, just about 0.3 - 0.5 MeV above the ground state energy of 
majority of configurations. Configurations corresponding to these local minima 
were characterized by a higher level of symmetry and possibly contained alpha¬ 
like cluster(s) in the middle. Surprisingly, binding energy distribution for these 
structures did show characteristic “kinks” around magic nuclei, despite some 
energy offset. These specific structures are those presented in Fig. 

The question that arose from the above observation was whether it might be 
possible to convert the local minima into global ones, while keeping the highly 
symmetric structures, by introducing the corrections to the Hamiltonian which 
take into account the cluster-cluster interactions. The following sections try to 
verify this hypothesis. 


4.1 Matter without alpha clusters 

The structures from Fig. |^have been used as structural constraints in the mini¬ 
mization procedure by providing the starting points for the preformed alpha-like 
clusters. In other words, the random initial configurations, as in Fig. have 
been replaced by those from Fig. The actual ground state configurations 
were determined through minimization of the liquid drop-like Hamiltonian ( |16[ ) 
by varying the mean positions and the widths of the individual wave packets. 

The following analysis is restricted to nuclei smaller or equal to ^°Ca which 
is the heaviest stable n-alpha nucleus. Beyond Z = 20 the nuclei become ra¬ 
dioactive and there are no data about their sizes available. 


The liquid drop-like Hamiltonian (161 with the parameters from Table 1 and 


Table 2 has been used together with structural constraints to obtain the ground 
state configurations for each of the considered nuclei. The minimization proce¬ 
dure provided the final parameters determining the ground state wave function: 
(r^) and (7^. Figure [^presents the obtained ground state binding energies and 
rms radii as a function of the atomic number Z. The experimental data and the 
model results are represented by the squares and the circles, respectively. As 
one can see the calculated binding energies change their slope at the same place 
{Z = 8) as the experimental ones and the same concerns the radii. The overall 
calculated binding energies are somewhat bigger by about 0.3 - 0.5 MeV/nucleon 
than the experimental ones and the calculated radii are somewhat smaller, how¬ 
ever both follow the experimental trends. Properties of particles up to "‘He are 
reproduced perfectly well, because they served as the constraints in the param¬ 
eter fitting procedure. 

In order to check the sensitivity of the results on the volume energy param¬ 
eter eoo similar calculations for eoo equal to -16.1 MeV instead of -16.79 MeV 
have been performed. The results are presented with a dotted green line in Fig. 
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Figure 5: The binding energies and rms radii of even-even nuclei (N = Z) ob¬ 
tained for the ground state wave functions and the Hamiltonian . The solid 
line connects the experimental points. Circles and triangles represent model 
results for two sets of parameters of the EoS. 


In case of the rms radii there is almost no difference observed. The binding 
energies still deviate from the experimental values, but what is important, sim¬ 
ilar changes of the slopes in theoretical and experimental trends are observed. 
Results of other trial calculations in which all the remaining parameters from 
Co to A were varied, did not yield better reproduction of the experimental val¬ 
ues with the use of the liquid drop form of the Hamiltonian together with the 
structural constraints i-iii. Better reproduction could only be obtained with the 
corrections to the Hamiltonian discussed in the next section. The structural con¬ 
straints assure configurations for which global minima of energy are obtained, 
provided that the alpha-alpha interactions are taken into account. Otherwise, 
they correspond to some local minima. 
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4.2 Modification of the Hamiltonian due to alpha clusters. 


Let us consider now formation of the nucleus by randomly located nucleons 
remaining within the range of their mutual interaction. Due to the Hamiltonian 
(16) nucleons are forming groups in such a way that the energy of the system is 
minimized. At certain density of the matter alpha-like clusters can be formed 
as a result of the spin-isospin interactions of nucleon groups. As alpha clusters 
are formed additional extra interactions connected with multi body effects can 
appear having a strong influence on the final structure of the system. Fig. 
illustrates a ground state of such a system. One can see that inside a matter 
drop alpha clusters appear and interactions between them are changing the 
energy of the system. Below, corrections to the Hamiltonian arising from these 
extra interactions will be defined. 



Figure 6: Visualization of the additional interactions between alpha clusters due 
to the many body interactions. 

Assuming that alpha-like clusters appear in the structure of the nucleus, one 
can expect two kinds of minor adjustments of the Hamiltonian, which should 
describe the essential part of the a — a interaction. 

The first correction corresponds to the possible interaction caused by the 
appearance of specific many-body forces. Such many-body effect can be inter¬ 
preted as a modification of the type and numbers of exchange bosons in the 
interaction process. It appears possible that, in such an exchange the charge 
currents are strongly limited because they change the energy of the considered 
clusters. In this case, we can expect additional term in the interaction potential 
which gives the repulsive forces between groups of nucleons forming a cluster of 
alpha-like particles (see Fig. |^. 

Here we test this hotfix in a form given by the product of the harmonic 
oscillator potential Vaa idij) ^.nd probabilities Pa (i ), Pa (j) that given four 
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nucleons (groups named i or j) form the alpha clusters: 


($ I AiJi I $) = ^ (i) (j) (d,,) 


( 18 ) 


where potential Vaa (dij) is taken in the form: 

^aa {.dij') S: {dij dg) 


(19) 


with dij being the average distance between clusters i and j, and Pq, (i) describ¬ 
ing the probability that a group of four nucleons 4n{i) forms an alpha cluster i. 
We assume that Pa {i) can be expressed as: 



( 20 ) 


where cr 4 „(j) is the variance of the distribution describing the probability of 
finding a nucleon in a group 4n (i) and aao describes such a variance (experi¬ 
mentally measured) for alpha particles. In the above formulas k, dg and are 
free parameters. It is assumed that the potential is equal to 0 when dij > 4 fm 
with a shape close to the Van der Waals one. 

The second correction is related to the additional symmetry of the wave 
function, which can be associated with the bosonic character of an alpha cluster. 

In the present approach nucleons and their clusters are distinguishable ob¬ 
jects Q. For nucleons their energy connected with the fermionic motion is taken 
into account by applying the energy density functional. For alpha particles a 
relevant correction to energy should also appear and it is an empirical one. 

When calculating the minimum of the Hamiltonian by varying the parame¬ 
ters of the wave function, the position of the alpha like cluster and its size and 
spread are established. 

These quantities can be different for different clusters depending on the 
adopted structural constraint for a given nucleus. The exchange of clusters must 
cause an increase of the energy of the system as a minimum of the Hamiltonian 
defines definite spreads of clusters for their different positions. 

For the indistinguishable alpha particles due to the symmetrization with 
respect to the exchange of cluster positions there is no chance to relate cluster 
sizes to their positions. 

The second correction is an empirical one: 



2 


( 21 ) 


where Na is the number of alpha clusters and /r is a parameter. <Ti{r) is the 
variance of the nucleon radius and crg(r) is the average variance of nucleons 
forming an alpha particle. 

Thus, this correction to the Hamiltonian describes an increase of the energy 
caused by different sizes of clusters and it should vanish when the wave functions 
for all the clusters correspond to alpha particles in their ground states. 
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It seems that this correction could be related to the bosonic character of 
alpha particles. Both corrections are positive and in this way the binding energy 
decreases. Assuming that parameters of the main part of Hamiltonian (14) are 
established by properties of light nuclei d, t, ^He and a, the idea now is to define 
additional four parameters in order to reproduce binding energies and sizes of 
nuclei starting from ®Be till "^^Ca by imposing the structural constraints. 

The adopted values of these parameters are summarized in Table 3. 


d 

V 

do 

K, 

[MeV] 

[1] 

fm 

MeV/ frri^ 

1. 

9. 

3.77 

5.2 


Table 3: Estimated parameters of corrections (18), and (21) to the Hamiltonian 
(16) coming from interactions between alpha clusters. 


The results for the full Hamiltonian, including the corrections due to the 
alpha-alpha interactions are presented in Fig. The upper part of the figure 
shows the binding energies as a function of the atomic number Z. The bottom 
part presents the rms radii. 

Much better agreement between the model results and the experiment can 
be observed now, as compared to Fig. Some discrepancies can be still noticed 
especially for ®Be for which the distance between the alpha particles is pretty 
small and much smaller than do in the harmonic potential description. Again, 
the properties of particles up to "^He are reproduced perfectly well, because they 
served as the constraints to fix the EoS parameters. 

Figure [^presents a comparison of the model and experimental density pro¬ 
files for the selected nuclei. The full (red) lines represent the experimental 
profiles and the broken (blue) lines result from the model calculations. Experi¬ 
mental density profiles were taken from j^, however for ^°Ne and ^®Ar reliable 
experimental data were not available. 

Figure]^ shows the isodensity contour plots for the model ground state con¬ 
figurations of ®Be, and ^°Ne. Please note that apart from the ^®0, 

which has a tetrahedral shape, all the other nuclei have a planar geometry 
(cf. Fi g. [4| . Similar shape for has been obtained within the FMD 28 and 
AMD |29]~simulations. In particular, the ^°Ne case shows that the central a-like 
cluster is more diffused than the surface ones. More quantitatively, the standard 
deviations of the density distributions of surface alphas are of the order of 1.02 
fm, while for the core alpha they are of the order of 1.28 fm. 

Figure [^presents a distribution of distances between alpha clusters for the 
nuclei in question. 

The distances between clusters from the ground state configurations obtained 
with the liquid drop-like Hamiltonian without corrections are represented by 
open circles. The distances obtained from the full Hamiltonian calculations are 
represented by closed ones. As can be seen, in the latter case the distances are 
grouping around 3.5 fm. This can denote that these clusters behave more like 
rigid spheres. 
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Figure 7: Binding energies per nucleon and rms radii for the Hamiltonian with 
the corrections (181 and (211. The right scale of binding energies per nucleon 
in the top panel refers to the magnified view for nuclei from to ^°Ca. The 
solid (dashed) line connects the experimental (model) points. 


5 Summary and Conclusions 

A simple model of nuclei formed out of alpha-like clusters in which interactions 
lead to strong correlations between nucleons has been presented. The model 
interactions have been defined using a cubic approximation of the nuclear EoS. 
The proposed EoS contains additional terms which take into account the spin 
and isospin polarizations. These additional terms cause that nucleonic clusters, 
with an enlarged binding energy, appear in the ground state configurations. 
These clusters have full spin-isospin symmetry similarly as in and in . For 
u-alpha nuclei these clusters can be identified as alpha like structures. When 
these clusters are formed in regions of smaller density, i.e. at the nuclear surface, 
then their sizes measured as the variance of the density distribution become 
comparable to the alpha particle size. On the other hand formation of the 
clusters in the central part of nuclei leads to much bigger spreads of their density 
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Radius [fm] 


Figure 8: Experimental (solid, red) and model (dashed, blue) density profiles. 



-5 0 5 -5 0 5 

X [fm] X [fm] 


Figure 9: Shapes of the model nuclei. The isodensity lines correspond to 1, 15, 
50, 75 and 95% of the maximum density. 
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Figure 10: The distribution of the distances between alpha clusters as a function 
of Z. The open (closed) circles represent distances without (with) the correction 
for alpha-alpha interaction in the Hamiltonian. 


distributions and they can hardly be identified as alpha particles. 

It has been shown that by taking into account the liquid drop-like form of 
the Hamiltonian in the minimization process it was not possible to describe the 
subtle changes of the binding energies of nuclei. Therefore the liquid drop part of 
the Hamiltonian has been supplemented by the corrections taking into account 
the cluster-cluster interactions. The ground state configurations obtained from 
the minimization of the full Hamiltonian with the imposed structural constraints 
were found to reproduce very well the experimental binding energies and sizes 
of the n-alpha nuclei and the available density profiles. In particular, the model 
with the cluster corrections has been able to reproduce the "fine structure" of 
the ground state energy distribution (inset of Fig. |^. In this respect it seems 
to be able to mimic the shell effects without explicit spin-orbit term. 

The interactions have been defined using a 12 parameter EoS approximation. 
The thirteenth parameter has been added to control the variance of the field 
for finite systems. The values of these parameters and correlations among them 
have been constrained using the experimental binding energies and sizes of d, 
f, ^He and a particles. Additional 4 parameters have been introduced in the 
corrections to the liquid drop part of the Hamiltonian in order to account for 
the interactions between the alpha-like clusters and to substantially improve the 
model predictions. The differences between the calculated and measured ground 
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state energies for nuclei with Z>8 are about an order of magnitude smaller than 
those obtained using the liquid drop model. Besides, the obtained parameters 
describe the properties of the t, ^He and a very accurately (practically error 
free). 

It has also been demonstrated that correlations between nucleons in the 
phase space can play an important role in the description of the ground states, 
particularly of the ground states of nuclei with the same and even number of 
neutrons and protons. 
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Appendix: Extraction of the EoS parameters 

The Hamiltonian (16) depends on 13 parameters and since the binding energies 
and rms radii of the ground states of d, t, ^He and a provide only 8 constraints, 
5 of the parameters have to be fixed by some other means. 

The first step of the procedure to fix the EoS parameters consisted in con¬ 
straining the values of Kq, pq, ego and A by the binding energy and rms radius 
of an alpha particle. This step has been summarized in the flow chart of Fig. 

HB 

The compressibility parameter of symmetric matter Kg and the volume en¬ 
ergy parameter eoo have been fixed within the range of the recent experimental 
constraints. The values of the remaining two parameters, pg and A are then 
uniquely determined from the correlation plot of Fig. The remaining 9 pa¬ 
rameters were irrelevant for the properties of an alpha particle, which greatly 
simplified the corresponding EoS. 

In the second step, the parameters e/o, Li and AT/ have been found in an 
analogous way by using the binding energies and the rms radii of t and ^He as 
constraints in eq. (17), by fixing the previously determined parameters Kg, pg, 
egg and A, and by disregarding the remaining 6 parameters. 

Finally, the remaining 6 parameters Bag, La, Ka, Cijg, Lij and K^j have 
been found by minimizing (17) with the fixed 7 parameters determined in the 
previous two steps and using the binding energies and the rms radii of d, t and 
^He as constraints. Fitting is done in an analogous way as in the inner loop of 


23 



Figure 11: Flow chart of the procedure of fixing the K^, po, eoo and A parameters 
using the binding energy and the rms radius of an alpha particle as a constraint. 
The details can be found in Sect. 3. 

the flow chart of Fig. El replacing the eoo and A parameters by the remaining 
6 ones. 

After fixing the 13 parameters of the liquid-drop-like Hamiltonian (16), the 
additional four parameters of the a — a interaction, p, v, do and k, have been 
fixed by fitting the ground state binding energies and rms radii of the nuclei 
from ®Be till "^^Ca. Here, in addition the structural constraints i. - iii. have 
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been imposed in the fitting routine, meaning that only the widths of the nucleon 
wave packets and the inter-a-like clusters distances could be varied, preserving 
the imposed configuration. 

The procedure of finding the ground state configurations which minimize the 
model Hamiltonian is described in Sect. 4 of [^, see also Fig. 
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